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GENERALIZATION OF THE FRACTIONAL POISSON 

DISTRIBUTION 

Richard Herrmann ^ 


Abstract 

A generalization of the Poisson distribution based on the generalized 
Mittag-LefHer function is proposed and the raw moments are cal¬ 

culated algebraically in terms of Bell polynomials. It is demonstrated, that 
the proposed distribution function contains the standard fractional Pois¬ 
son distribution as a subset. A possible interpretation of the additional 
parameter (3 is suggested. 
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1. Introduction 

Since its early beginnings, fractional calculus has developed as an ever 
growing powerful tool to model complex phenomena in different branches 
of e.g. physics or biology. 

A major field of research is the area of fractional Brownian motion and 
anomalous diffusion phenomena. Closely related are random walk problems 
and the investigation of stochastic processes, where e.g. survival probabili¬ 
ties R{t) differ from the exponential law, which is known from the classical 
Poisson process. 
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Based on the concept of non-exponential power law Markov-processes, 
the fractional extension of the classical Poisson process leads to the stan¬ 
dard fractional Poisson distribution function of the form I8],[I2],[9],[I3],[6] 


Pk{a,X) 


kl ^ re! r(a(/c-|-re)-|-1) 

n=0 

spit, .11-,.. 

(-1)‘(^^£«(A))I-„a. 0<a<l (1.3) 


which may be considered as a generalization of the classical Poisson distri¬ 
bution represented by the introduction of the Mittag-Leffler function Ea{X) 
and its derivatives |10j . 

The mean /r and the variance a for this fractional Poisson distribution 
follow as [6] 


a 




r(i + a) 


2/ 0rr(l-ha) 


P + P 


22a-lr(i/2 + Q,) 


-1) 


(1.4) 

(1.5) 


The distribution function dnD is the result of the underlying process 
and allows for a direct comparison of experimental data with the corre¬ 
sponding model assumptions. 

In the following, we will consider the inverse problem: We will present 
fractional generalizations of distribution functions, based on generalized 
Mittag-Leffler functions. 

The motivation for this approach is straight forward: A given exper¬ 
imental setup will generate distribution data first, which may be directly 
compared with generalized distribution functions. The underlying micro¬ 
scopic process may be investigated in a second step. 

A large amount of literature covers extensions of the standard Mittag- 
Leffler function by introducing additional parameters and functional be¬ 
havior (see e.g. [5], HI and references therein), but a direct e.g. physical 
interpretation of the additional parameters is a still open question in many 
cases. 

In the following we therefore will apply an approach, which allows for 
a wide class of Mittag-Leffler type functions to determine the correspond¬ 
ing distribution functions. This distributions may be directly compared 
with observed phenomena, which will turn out to be helpful for a broader 
understanding of the functionality and scope of possible applications. 
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2. Nomenclature 

In order to calculate the major properties of a generalized Poisson dis¬ 
tribution based on the generalized Mittag-Leffler function, let us first recall 
the major properties of the classical Poisson distribution, which is given 
by: 

«(A) = (2.1) 

with the normalization constant N, which is determined by the requirement 
of normalizability of the distribution, which coincides with the zeroth raw 
moment hq: 

^ 1 ^ 1 

k=0 k=0 

Higher raw moments are then given as: 

OO CXD . ^ 

Mn = A;fc>fc(A) = e-^ = (2.3) 

k=0 k=0 

where (12.3|] is the dehning equation for the Bell polynomials B^, which 
obey the recursion relation [T] 

-^Bnix) = -Bn+l{x) - Bnix) (2.4) 

ax X 

and are related to the Stirling numbers of second kind S{n,k) [T6| via 
Dobrinski’s formula[2]: 

n 

Bn{X) = Y.S{n,k)X^ (2.5) 

k=l 

The hrst raw moments result as: 

IB = X ( 2 . 6 ) 

M2 = A(A -|- 1) (2.7) 

and therefore the mean ^ and variance a follow as centered moments: 

= //I = A (2.8) 

a = )U2 - Ml = A(A -b 1) - A^ = A (2.9) 

3. Generalized fractional Poisson distribution 

We now extend the definition of the classical Poisson distribution pk{X) 
by applying the canonical fractionalization procedure: 
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We replace the faculty in (12.11) by the T function, introduce as an ex¬ 
ample two additional parameters a, {3 and define a generalized fractional 
Poisson distribution pk{X,a,f3) via; 


Pk{X,a,l3) 


1 

Nr{f3 + ak) 


(3.1) 


the normalization constant N follows from the requirement of normaliz- 
ability of the distribution: 


oo ^ oo 

Po = ^Pfc(A,a,^) = 




1 


k=0 


N + ak) N 




(3.2) 


with the generalized Mittag-Leffler function Ea^p{X) which is defined as 

A^ 


Ea.fW = Y. 


k=0 


r(/3 -|- ak) 


The distribution is therefore given by: 

Pk{X,a,(3) = 


1 


A^ 


E^^p{X)T{(3 + ak) 


(3.3) 


(3.4) 


Higher raw moments may be calculated by an iterative procedure, which 
we sketch for the case n = 1. 

Starting with unnormalized zeroth moment, we apply the operator ^A: 


EaA>^) = £ 


A^ 


k=0 

OO 


r(/3 -|- ak) 


-X-E^A^) = + 


A' 


Using the commutation relation 


fc=o 


r(/3 -|- ak) 


r d , , 

the result for the first raw moment follows as: 

/-7 OO » ^ 

X—■ E^A>^) + E^A>^) = 


+ E 


6.0 EW + ak) 6^r{/3 + Q.t) 


= E'‘ 


A‘ 


fc =0 


r(/3 -|- ak) 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 


This procedure may be applied iteratively n times to obtain the n-th raw 
moment. 
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We obtain an iterative solution for the raw moments: 


1 


— 


-(A 


d 




(3.10) 


dX' 

which results in a n-fold application of the Euler-operator Je = 

In order to deduce a closed form solution, we perform a normal ordering 
of the n-fold Euler operator 
With the settings 


X 


= {:X±:rE^,^iX) 


(3.11) 

(3.12) 


where E'^"^(A) denotes the n-th derivative of the Mittag-Leffler function 
with respect to A and (::) indicates the normal ordered product m 


(:A 


d 


= A^ 


d" 


dX ' dX^ 

the iteration scheme is now isomorphic to the recurrence relation 
we obtain for the raw moments 


— 




(3.13) 
and 


(3.14) 


En 


- k=n 

\2Sin,k){X^E^^],{X)) (3.15) 

^ k=i 


The mean /r and variance a follow as the lowest centered moments: 

^ ' - 7 ( 1 ) 


n = m = 


a = ^i2- lii = 


A^log(E„,^(A)) 


1 


(a2e(2)(a) + ae^;;,(a))-^^ 


(1) 


d 


(3.16) 

(3.17) 

(3.18) 


(A|-log(Ei”{A)))(ALlog(E„_„(A))) 

-(aL log(E„,^(A)))" + ALl„g(£:„ „(A)) (3.19) 


Optionally, with 

d 


^^E'Q, ,g(x) — ^ Eq, (x) -|- 


1-/3 


a 


Ea,a+f} ( 2 :) 


(3.20) 


all calculations may be reduced to Mittag-Leffler functions directly. 
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Figure 1. For A = 5 and /? = 1 the probability distribution 
Pfc(A = 5,a,/3 = 1) from (13.41) is plotted for 1.5 > a > 0.5 
in 0.05 steps. Thick dashed line indicates a = 1, which 
corresponds to the classical Poisson distribution. 


4. Discussion 


In contrast to the standard fractional Poisson distribution, which is 
valid only in the interval 0 < a < 1, our derivation is not restricted to a 
specific region of allowed a values. Nevertheless, in the following we present 
results in the practically interesting region 0 < a < 2. In figures [Hand [2] we 
have plotted the generalized Poisson distribution (13.411 for varying a and /3 
respectively. 

To simplify an interpretation, we first calculate the mean and variance 
for two limiting cases: 

First, for small A <C 1 we obtain 


lim u 
A^o 


lim a 
A^O 


nP) y 

T{a + P) 

r(/?) 

r(a + /3) 


(4.1) 

(4.2) 


which shows two interesting features: 
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k 

Figure 2. For A = 5 and a = 1 the probability distribution 
Pfc(A = 5, a = 1, /3) from (j3.4jl is plotted for 4 > /3 > —4 in 
0.4 steps. Thick dashed line indicates /? = 1, which corre¬ 
sponds to the classical Poisson distribution. 


The feature /r = u, which we know from the classical Poisson distribu¬ 
tion, holds for small A for the generalized case too. For /3 = 1 the mean 
corresponds to the mean value, which is obtained with standard fractional 
Poisson distribution only the shape differs (since variances are differ¬ 

ent). 

However, this case (A <C 1) is a very special case, in general, we deal 
with A ^ 1. 

For the case A S> 1, using the asymptotic expression for the generalized 
Mittag-Leffler function [3] 

lim F„fl(A) = +o(|) 0<a<2 (4.3) 

A^oo a A 


we obtain 


lim fi 

\—>-oo 


1-(3 + AV^ 
a 

AV« 


(4.4) 


lim a 
A—^oo 


(4.5) 
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standard versus 
generalized 
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Figure 3. Fit results of the standard fractional Poisson dis¬ 
tribution pkiois, -^) H-ip with u = 1 for 0.1 < Os < 1-0 in 0.1 
steps (dots) with generalized fractional Poisson distribution 
from pk{a,l3,X) (13.41) (lines) for A = 5. Fit parameter sets 
{a./3} are listed in table [Ij 

which allows for a direct interpretation of parameter /3: For a given a, the 
expectation value p is shifted by the amount (1 — (3)1 a, while the variance 
is independent of j3. In figure El this behavior is easily observed. 

Consequently, the parameter f3 could be directly interpreted within 
predator-prey systems: Since the probability distribution is shifted by a 
(3 depended constant amount from the expected mean ~ ja, while 
other characteristics (e.g. shape) of the distribution remain unchanged, /3 
may be considered as a constant time independent population decrease and 
increase rate respectively (comparable to a sink/source in electro dynam¬ 
ics). 

With this understanding, in a next step we may then investigate the 
corresponding fractional process, considering a fractional pendant of the 
coupled sets of Lotka-Volterra equations [H], which are currently used to 
model such scenarios. 









GENERALIZED ERACTIONAL POISSON DISTRIBUTION 9 


Table 1. Parameter sets for a fit of the standard frac¬ 
tional Poisson distribution pf:{as,X) (jl.ip with z/ = 1 for 
different Ug and the generalized fractional Poisson distribu¬ 
tion Pkia, /?, A) (|3.4p for A = 5. 


as 

a 

/? 

(Xg 

a 

/3 

1.0 

1.000 

1.00 

0.5 

0.755 

10.40 

0.9 

0.911 

2.86 

0.4 

0.718 

13.73 

0.8 

0.855 

4.58 

0.3 

0.670 

19.40 

0.7 

0.815 

6.33 

0.2 

0.606 

31.61 

0.6 

0.781 

8.26 

0.1 

0.507 

76.01 


Einally, the parameter set {a, j3} in (13.411 may be adjusted appropri¬ 
ately to obtain similar mean p and variance a values of both, the standard 
fractional Poisson distribution (ini) and the generalized fractional Poisson 
distribution , which leads to similar shapes of the distribution function. 
This may be realized by calculating the parameter pair {a, /?} by equating 
p and a from (11.41) . (11.51) with (13.161) . (13.181) . by the requirement, that both 
probability distributions should be same values for two given k’s or by a fit 
procedure, e.g. near the maximum of the distribution. 

In figure [3] we compare the results of a fitting procedure of the fractional 
Poisson distribution (jl.ip with (j3.4p . In tabled] the corresponding fitted pa¬ 
rameters are listed. The agreement is very good. In practical applications, 
it would be difficult, to distinguish both distribution families. 

Therefore we may conclude, that the proposed generalized Poisson dis¬ 
tribution covers the complete parameter range of the standard fractional 
Poisson distribution as a subset. It should at least be mentioned, that the 
numerical behavior of the proposed generalized fractional Poisson distribu¬ 
tion is much easier to handle than the standard distribution. 


Conclusion 

We have presented a general method to derive the characteristic prop¬ 
erties of a given fractional distribution in terms of normal ordered Bell 
polynomials with the Euler-operator as an argument. 

We have demonstrated for the special case of the generalized Mittag- 
Leffler function, that the investigation of corresponding fractional distri¬ 
butions helps to obtain a clearer understanding of first the meaning of 
additional parameters in generalizations of Mittag-Leffler type functions 
and second of possible areas of their application. 
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In addition we have shown, that the standard fractional Poisson distri¬ 
bution may be reproduced as a specific subset of the generalized fractional 
Poisson distribution. 
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